Cosmological Constraints from 21cm Surveys After Reionization 
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21cm emission from residual neutral hydrogen after the epoch of reionization can be used to trace 
the cosmological power spectrum of density fluctuations. Using a Fisher matrix formulation, we 
provide a detailed forecast of the constraints on cosmological parameters that are achievable with this 
probe. We consider two designs: a scaled-up version of the MWA observatory as well as a Fast Fourier 
Transform Telescope. We find that 21cm observations dedicated to post-reionization redshifts may 
yield significantly better constraints than next generation Cosmic Microwave Background (CMB) 
experiments. We find the constraints on Qa, Qmh 2 , and Q. v h 2 to be the strongest, each improved by 
at least an order of magnitude over the Planck CMB satellite alone for both designs. Our results do 
not depend as strongly on uncertainties in the astrophysics associated with the ionization of hydrogen 
as similar 21cm surveys during the epoch of reionization. However, we find that modulation of the 
21cm power spectrum from the ionizing background could potentially degrade constraints on the 
spectral index of the primordial power spectrum and its running by more than an order of magnitude. 
Our results also depend strongly on the maximum wavenumber of the power spectrum which can 
be used due to non-linearities. 

PACS numbers: 



I. INTRODUCTION 

Recently, there has been much interest in the feasi- 
bility of mapping the three-dimensional distribution of 
cosmic hydrogen through its spin-flip transition at a res- 
onant rest frame wavelength of 21cm [H, 0|. Several first 
generation experiments are being constructed to probe 
the epoch of reionization (MWA [|, LOFAR [|, PAPER 
[|, 21 CM A 0) and more ambitious designs are being 
planned (SKA 0). 

One driver for mapping hydrogen through 21cm emis- 
sion is to measure cosmological parameters from the un- 
derlying cosmic power spectrum. During the epoch of 
reionization (EoR), the 21cm power spectrum is shaped 
mainly by the structure of ionized regions. Even without 
precise knowledge of the ionization power spectrum it is 
possible to isolate the cosmological power spectrum by 
exploiting anisotropics in redshift space due to peculiar 
velocities @,@|. 

Recent work [13, EH has shown that the 21cm power 
spectrum accessible during the EoR has the potential to 
put tight constraints on cosmological parameters; how- 
ever, these constraints depend on model-dependent un- 
certainties, the most important of which involves the ion- 
ization power spectrum. 

In this paper, we explore in detail the constraints 
achievable from map ping the residual cosmic hydrogen 
after reionization [ll EE El EE EE E3 • The hydrogen 
resides in pockets of dense galactic regions which are self- 
shielded from the UV background, also known as damped 
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Lyman-cv absorbers (DLAs) in quasar spectra [l8| . In 
difference from traditional galaxy redshift surveys, 21cm 
surveys do not need to resolve individual galaxies but 
rather are able to monitor the smooth variation in their 
cumulative 21cm emission owing to their clustering on 
large scales fl5l |. 

Measuring the 21cm power spectrum after reioniza- 
tion offers several key advantages as compared to the 
EoR. First, the near uniformity of the UV radiation field 
guarantees that the 21cm power spectrum would reliably 
trace the underlying matter power spectrum. However, 
the UV background could introduce a scale dependent 
modulation to the 21cm power spectrum as large as one 
percent [ijj • This modulation could introduce degenera- 
cies that would significantly weaken constraints on the 
spectral index of the primordial power spectrum and its 
running. Another advantage of probing low redshifts 
is that the brightness temperature of the galactic syn- 
chrotron emission scales as (1 + z) 2 ^ . However, this ad- 
vantage is offset by the fact that the mass weighted neu- 
tral fraction of hydrogen is only a few percent at redshifts 
z < 6 PJ]. 

In this paper we use the Fisher matrix formalism to 
quantify how effectively futuristic surveys dedicated to 
post-reionization redshifts (z < 6) can constrain cos- 
mology. We consider both a survey with ten times the 
number of cross dipoles of the MWA 0] (5000 x 16 = 
8 x 10 4 ) but operating at higher frequencies, which we 
term MWA5000, and a Fast Fourier Transform Telescope 
(FFTT) with 10 6 dipoles over a square kilometer area 
[201 ] . We also show results for the combination of these 
surveys and a next generation Cosmic Microwave Back- 
ground (CMB) experiment (Planck). 

This paper is structured as follows. In §11 we discuss 
the detectability of the 21cm signal after reionization and 
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in §111 we discuss its power spectrum. In §IV we consider 
the details of our Fisher matrix calculation and in §V we 
present its results. Finally, we discuss and summarize 
our conclusions in §VI and §VII. 



to these redshifts to obtain the type of results discussed 
in this paper. 

III. 21CM POWER SPECTRUM 



II. DETECTABILITY OF THE 21CM SIGNAL 
AFTER REIONIZATION 

Reionization starts with ionized (HII) regions around 
galaxies which grow and eventually overlap. This overlap 
of HII regions characterizes the end of the EoR. It was 
once thought that the 21cm signal would disappear after 
this transition because there is little neutral hydrogen 
left. Recent work has shown that this is not the case 
CDI [H HI E2|. The detectability of the signal may 
not decline substantially following the end of the EoR 
because the Galactic synchrotron foreground is weaker at 
higher frequencies (lower redshifts of the 21cm emission). 

After reionization most of the remaining neutral hy- 
drogen is expected to reside in DLAs. Observations have 
shown that out to z w 4 the cosmological density pa- 
rameter of HI is fiHi ~ 10 -3 [19]. This corresponds to a 
mass-averaged neutral fraction of a few percent. It does 
not contribute significantly to the Lya forest, which is 
mostly shaped by the much smaller volume-averaged neu- 
tral fraction. A 21cm survey like the ones we discuss are 
sensitive to all neutral hydrogen within the survey vol- 
ume because no galaxies are individually identified and 
thus there is no minimum threshold for the detection of 
individual galaxies. 

Even though the majority of the neutral hydrogen re- 
sides in self shielded clumps, self absorption is not ex- 
pected to significantly reduce the 21cm signal. This is 
supported by 21cm absorption studies of DLAs over the 
redshift interval < z < 3.4, which exhibit an optical 
depth to 21cm absorption of the radio flux from a back- 
ground quasar of less than a few percent [23|, [24[ ■ These 
observations are supported by theoretical calculations of 
the 21cm optical depth of neutral gas in high redshift 
minihalos [25|]. We also note that the spin temperature 
in DLAs is much higher than the CMB temperature and 
thus the 21cm signal is independent of the kinetic tem- 
perature of the gas [H[ . 

At z w 4 the 21cm brightness temperature contrast 
with the CMB will be roughly 0.5 mK. On scales of 10 
comoving Mpc, the rms amplitude of density fluctuations 
is roughly a w 0.2, so we expect 21cm fluctuations at a 
level of ~ O.lmK. This is only an order of magnitude or 
so less than the largest fluctuations expected during the 
EoR [26|. These fluctuations combined with the lower 
brightness temperature of galactic synchrotron emission, 
which scales as (1 + z) 2 ' 6 , should provide a detectable sig- 
nal. Previous work (e.g. Fig. 10 in Ref.(22j) shows that 
the signal to noise ratio of the 21cm signal after reioniza- 
tion may be similar or even higher than that during the 
EoR. However, it is important to note that it will be nec- 
essary to build special instruments which are optimized 



Next, we describe the 21cm power spectrum after the 
epoch of reionization. A complete discussion of the rel- 
evant physics can be found in Refs. [J, [H, [TH, [It}- The 
difference between the average brightness temperature of 
21cm emission at redshift z and the CMB temperature 
is described by 
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where xn is the global neutral hydrogen fraction and T s 
is the HI spin temperature. In Fourier space, the power 
spectrum of 21cm brightness fluctuations is defined by 
(27r) 3 (5 3 (k-k')P2i(k) = (ATb(k)*AT b (k')>, where k is 
the wave- vector of a given Fourier mode and ATb(k) is 
the brightness temperature fluctuation in Fourier space. 

Before giving its exact form, we stress that the most 
important feature of f2i(k) after reionization is that it 
essentially traces the cosmological matter power spec- 
trum, P{k). Extracting the cosmological power spectrum 
will make it possible to place tight constraints on many 
cosmological parameters. During the EoR the situation 
is more complicated. P2i(k) contains additional terms 
which depend on the ionization power spectrum and the 
cross power-spectrum between the ionization and matter 
distributions. 

Following the derivation of Ref. [l~7| we find that after 
reionization 



P 21 (k)=T 2 x 2 m P(k) [S(k) + //i 2 



(2) 



Here f h = T h /(x m (T s - Tcmb)/T s ), B(k) is the average 
scale dependent, xm mass- weighted bias for the DLA 
host galaxies, and /i = cos 9 where 9 is the angle between 
the line of sight and the wave- vector k [l7l . ui\ . The fac- 
tor [B (k) + f /j, 2 ] 2 arises due to redshift space anisotropics 
from line of sight peculiar velocities 8, 27, 28]. Here the 
growth index is defined as / = d\nD/d\na, where D(z) 
is the growth factor of density perturbations and a is 
the cosmic scale factor. The bias factor arises because 
the residual neutral hydrogen is located in host galax- 
ies which are biased with respect to the underlying dark 
matter distribution. This bias and the neutral hydrogen 
fraction are the only quantities in the above equation 
which do not depend solely on fundamental physics. 

While the near uniformity of the ionizing UV back- 
ground causes P2i(k) to trace the cosmological power 
spectrum, fluctuations introduce a modulation that could 
modify the power spectrum by less than one percent [I?) • 
This modulation of the power spectrum adds a scale de- 
pendent correction factor to the host galaxy bias. The 
average bias B(k) is equal to the mean halo bias 6, mul- 
tiplied by a factor (1 — K{k)). In appendix A, we derive 
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an expression for K(k) for the case where the mean-free- 
path of ionizing photons is independent of position. If 
the effect of the ionizing background is near the max- 
imum which could be expected it may be necessary to 
model it in order to estimate cosmological parameters 
reliably. In §V, we show how cosmological constraints 
would likely be affected by degeneracies with nuisance 
parameters introduced by a scale dependent bias model. 

The bias, b, for halos of a particular mass M can be 
derived using the Press-Schechtcr formalism modified to 
include non-spherical collapse [2!| 



was due to reionization modifying the average stellar age 
of galaxies and hence the mass to light ratio. We do 
not expect this to significantly effect the neutral hydro- 
gen hosted in galaxies. Additionally the modulation of 
the power spectrum was found to be approximately scale 
independent, which would not be degenerate with cos- 
mological parameters discussed in this work. 



Non-Linear Effects on Baryon Acoustic 
Oscillations 



b{M, z) = 1 
1 



+ bv' 2 ^- 



,/2c 



/'2c 



+ 6(l-c)(l-c/2) 



(3) 



where v' 2 = aS 2 /a 2 (M, z), a = 0.707, b = 0.5, and c = 
0.6. Here a 2 (M,z) is the variance of the density field 
smoothed on a mass scale M and S c is the linear over- 
density threshold for collapse at redshift z. 

Assuming that the neutral gas to halo mass ratio is 
independent of mass, we can derive the mean bias at a 
particular redshift, 6, as the weighted average 



M„ 



—MdM 



(4) 



where dn/dM is the Sheth-Tormen mass function of dark 
matter halos [291 ]. The minimum halo mass M m i n de- 
fines the threshold for assembling heated gas out of the 
photo-ionized intergalactic medium, corresponding to a 
minimum virial temperature T V i r w 10 5 K (30. [3ll. [32. l33l. 
H3, HH . The appropriate value of M m ; n will depend on 
the mass of galaxies which host DLAs, however we find 
our results insensitive to the exact choice of M m - m (see 
Fig. ID). 

A constant virial temperature leads to the redshift de- 
pendence M min = 5.3 x 10 9 ((1 + z)/4.5) -15 M Q . Based 
on observations of DLAs [l9| , we assume a fiducial value 
of neutral hydrogen density to be cchii = 0.02. 

Previous work has suggested that at high redshifts the 
gala xy power spectrum may be affected by reionization 
|36l Wf\ . In [36| it is shown that fluctuations in the red- 
shift of reionization in different locations leads to varia- 
tion in the minimum mass of galaxies which can be as- 
sembled. This creates a scale dependent modulation of 
the power spectrum. Recent observations find that at z ss 
3 galaxies with virial velocities < 50 — 70km/s contribute 
little to DLAs and that most DLAs are hosted within ha- 
los having viral velocities between V = 100 — 150 km/s 
[38| . Since DLAs and thus most of the remaining neutral 
hydrogen after reionization is not found in the smallest 
galaxies we do not expect changes in the minimum mass 
of galaxies which are assembled to have an important 
effect on the 21cm power spectrum. 

In [37[ , the effects of reionization on the galaxy power 
spectrum for more massive galaxies are explored. Reion- 
ization is found to change the power spectrum, but this 



It has been recently shown that non-linear effects may 
change the baryon acoustic oscillations (BAO) signature 
at larger scales than those previously considered [39| . 
The erasure of BAO information due to non-linear effects 
can be written as 



fb,lin(fe)x 
k 2 



exp 



^((l-A^Ei + ^Ej) 



(5) 



where P\> is the part of the power spectrum which con- 
tains the wiggles from the BAO, Ej_ = T,qD(z), Eii = 
S (l + f)D(z), E = ll^/i^Mpc, and D(z) is the 
growth function normalized to (1 + z) -1 at high redshifts. 
In our Fisher matrix calculations (Pb,ni(fc) + -Psmooth(fc)) 
is used in place of P(k) in Eq. [3J Pb.iin(fc) = P(k) — 
-Psmooth(fc), where P smoo th(&) is a low order polynomial 
fit to P(k) in log-log space. This fit removes the wig- 
gles from the BAO, such that P^ is just the part of the 
power spectrum containing the BAO wiggles. Our choice 
of E corresponds to the value determined by simulations 
and scaled linearly with erg to fit our fiducial cosmology. 
The scales implied by E^ and Em are held fixed in ob- 
servable units such that the exponential suppression can 
be taken outside of the derivatives in the Fisher matrix 
formalism described below. This ensures that the BAO 
erasure does not in any way improve the forecasted errors 
on cosmological parameters. 

We emphasize that the exponential suppression of the 
power spectrum in Eq. [5] is only applied to the the 
BAO wiggles after the smooth linear power has been 
subtracted off. The power spectrum used in the deriva- 
tives which enter the Fisher matrix described below is 
the sum of the wiggles suppressed by Eq. [5] and the un- 
suppressed smooth linear power spectrum. For the error 
on P21 which enters the Fisher matrix we do not sup- 
press either the wiggles or the smooth part of the power 
spectrum. 



IV. FISHER MATRIX FORMULATION 

Given a set of parameters Ai, the Fisher matrix for- 
malism provides an estimation of the error for each of 
the parameters associated with some data set 40] . The 
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la errors on parameters can be estimated as 



AA, : 



where 



Fij — 



E 

pixels 



l 



SP 2 \ d\ 



dp 



21 



dP, 



21 



d\ 4 



(6) 



(7) 



9, 10], P21 is the total 21cm power spectrum, and SP21 is 
the uncertainty on a measurement of the power spectrum. 
The relevant derivatives have been calculated using the 
transfer functions from CAMB [411 ]. 

We work in u-space rather than k-space as described 
in Refs. [ij [l(| to simplify calculating Alcock-Paczynski 
effects. A radio interferometer directly measures visibil- 
ities, 



V{u,v,v) = / dnAT b (n,j/)A„(n)e 



(8) 



where V is the visibility for a pair of antennae and A v is 
the contribution to the primary beam in the n direction. 
Here we have used the flat sky approximation. This is ap- 
propriate even in the case of the FFTT, which images the 
entire hemisphere, because essentially all of the cosmolog- 
ical information in our surveys is found on small angular 
scales. The vector = (u, v) corresponds to the number 
of wavelengths between the pair of antennae. Performing 
the Fourier transform I(u) — J dvV{u, v, v) exp(27rii/?7), 
we obtain a signal in terms of u = ui + vj + 77k, where 
77 has units of time and k is the unit vector along the 
line of sight. Note that there is a one to one correspon- 
dence between k and u given by 2ttu±/<1a = k^ per- 
pendicular to the line of sight and 27rui|/y = k|| along 
the line of sight, where d,A is the angular diameter dis- 
tance to the observation and y is the ratio of comoving 
distance to frequency interval. With the Fourier conven- 
tions above we have p2i(k) = P2i(u)d A y. Working in 
u-space simplifies our calculation because P2i(u) is mea- 
surable without cosmological assumptions and thus Eq. 
(|7|) can be applied directly. If we were to work in k- 
space, the Alcock-Paczynski effect would distort P2i(k) 
from Eq. |2]| when we take derivatives with respect to 
cosmological parameters. 

Because we are considering discrete sources, the total 
21cm power spectrum will be the sum of Eq. |(5J) and 
a shot noise term P s hot = P2i(k)/(b 2 P(k)riY)LA), where 
"dla is the effective number density of the galaxies which 
host DLAs, P2i(fc) is P2i(k) along the line of sight, and 
b 2 P(k) is the power spectrum of DLA hosting galaxies 
along the line of sight. 

Any constant |k| and 9 define an annulus of constant 
p2i(u). The Fisher matrix is calculated from Eq. ([7|) by 
summing the contribution from annuli which fill all of 
k-space accessible to an observation. We divide a 21cm 
survey into redshift bins small enough for the redshift 
evolution of P(k) across a bin to be negligible. A sep- 
arate Fisher matrix is calculated for each redshift bin 



and then summed into a Fisher matrix which reflects the 
information about the entire survey. 

To calculate the error in Eq. ([7]) we follow the work of 
Refs. Hill. 

The error on P21 for a particular |k| and 9 

is 



SP21 



P21 + P 



shot 



P 



N 



y/N c x AW x P to t/P' 



(9) 



where P/v is the noise power spectrum of an interferom- 
eter, B is the bandwidth over which foregrounds can be 
removed, B tot is the total accessible bandwidth, N^ c \^ s is 
the number of fields being imaged, and N c is the number 
of independent cells in k-space. The number of cells in 
an annulus of constant P 2 i(u) is well approximated by 



AT r = 2irk 2 sin9AkA9- 



V 



(27T) 



(10) 



where V is the comoving volume observed by the exper- 
iment, A9 and Ak are the angular and k-space widths 
of the annulus. We choose Afc and A9 such that P21 (u) 
is essentially constant. The noise of an interferometer is 
given by 



Pn(ux) = 



A 2 P 



sys 



1 



t n(u±) ' 



(11) 



21cm x (1 + z) is the observed wavelength, 



L Sys is the system temperature of the interferometer, A c 



where A 
P 

is the effective area, rj is the total observing time, and 
n(u±) is the number density of baselines. This error is 
identical to that applicable in studies of the EoR with 
the exception of the shot noise [1, [l(| ■ 

The sensitivity of 21cm measurements will be impacted 
by how effectively foregrounds such as the galactic syn- 
chrotron emission can be removed. Refs. [§, H^, 44 1 
suggest that after fitting out a low order polynomial, 
residual foregrounds in the power spectrum will be neg- 
ligible if the frequency band over which foregrounds are 
removed, P, is substantially smaller than the total band- 
pass available. This is true for wave-vectors greater than 
&min = 27r/ (yP). Recent work estimating errors on cos- 
mological parameters with 21cm surveys during the EoR 
has used this assumption in modeling foreground removal 
0, EH El • It is encouraging that we find parameter con- 
straints to depend relatively weakly on the exact value of 
the bandwidth over which foregrounds are removed (see 
Fig. [2]) . We note that while we expect these assumptions 
pertaining to foreground removal to be robust, the true 
ability to remove foregrounds may not be known until 
first generation 21cm instruments are operational. 

There are several factors which limit the scales acces- 
sible to 21cm surveys for parameter determination. As 
discussed above, we assume that foregrounds can not be 
effectively removed from pixels with k < fc m j n . There is 
also a minimum accessible value of k sin 9 imposed by the 
minimum baseline, 6 m i n « \/ A p \ lys . We exclude informa- 
tion for k values above the linear scale. We compare 
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FIG. 1: Maximal wavenumber, fc max , above which the 
non-linear power-spectrum deviates from the linear power- 
spectrum by more than 25%, 20%, 15%, 10% and 5% (top to 
bottom). Larger discrepancies produce higher values of fc max . 



the linear power spectrum for our fiducial cosmology to 
a non-linear power spectrum produced by HALOFIT in 
CAMB [4lJ, and define our non-linear cutoff as the k 
value where there is a 10% discrepancy. 

The precision of the measurements discussed in this 
paper will require a theoretical prediction of the mat- 
ter power spectrum which is accurate to better than one 
percent. The k cutoff is meant to represent the scale to 
which future N-body simulations, such as those discussed 
in [ifl Ull , will be able to calculate the power spectrum 
to the required accuracy. Since the exact value of this 
scale is unknown we present the dependence of varying 
it on cosmological parameter constraints in Fig. [3l 

The value of this cutoff is fairly conservative when 
compared to the criterion of Ref. [47} which chooses 
kni = n/2R for a value of R giving an rms density fluc- 
tuation amplitude a(R) = 0.5. The latter criterion has 
been used for 21cm parameter constraints at higher red- 
shifts than we explore [TlT ] . The fc max values correspond- 
ing to different levels of discrepancy between the linear 
and non-linear power spectra are shown in Fig. [T] 

We also include constraints when combined with future 
data from the Planck satellite [4q . The method used to 
calculate the corresponding Fisher matrix was adopted 
from Refs. [Mill- 



A. Survey Characteristics 

We apply the Fisher matrix formulation to future 
surveys dedicated to low redshifts similar in design 
to MWA5000 and a Fast Fourier Transform Telescope 
(FFTT). We show results with surveys designed for a 
central rcdshift of 3.5 and 1.5 that each span a factor of 
3 in (1 + z). 

We model the MWA5000 observatory as an interferom- 



eter with 5000 tiles each containing 16 dipole antennae. 
For the survey centered at z = 3.5(1.5) we assume a con- 
stant core of antennae out to 40(22)m and then a r~ 2 
distribution out to 580(278)m. This layout gives roughly 
the same baseline density distribution as MWA does at 
higher redshifts [?|. 

The angular resolution of the telescope is set by the 
baseline density distribution. The longest baseline cor- 
responds to the highest value of /cj_ that can be probed, 
2iTVL± tmax I d,A = k±. max . The above distribution of an- 
tennae contains baselines which are long enough to re- 
solve k max at all redshifts. We assume a frequency res- 
olution of 0.01 MHz for all surveys considered. These 
arrays will have better resolution than this, but this is 
sufficient to resolve the scales set by k max . 

For the survey centered around z — 3.5 the effective 
area of each tile A e w 7Vdi p A 2 /4 and is limited by the 
physical area of the tile, A p h ys [5l|. The values of A p h ys 
are chosen such that Aphys ~ A c at the central redshift. 
For the low-redshift case, we assume that dishes are used 
instead of antennae, resulting in A c sa A p h ys . 

We model an FFTT observatory [2(| as an interferom- 
eter with 10 6 evenly spaced dipoles over a square kilo- 
meter that can all be correlated. The noise in each 
dipole is calculated separately with A c ss A 2 /4 limited 
by Aphys = S 2 , where S is the spacing of the dipoles on 
a square grid. 

We have assumed that foregrounds can be removed on 
scales of up to 30[4.5/(l + z)]MHz. This corresponds 
to the value that [15[ cite as desirable to measure the 
neutrino signal. The exact choice of this bandwidth is 
relatively unimportant, as our results depend weakly on 
this value (see Fig. [3]). We assume that 0.65 of the sky 
in one hemisphere is imaged for 2000 hours. This corre- 
sponds to roughly 16 fields of view. The partial coverage 
is due to foregrounds which can not be removed in the 
vicinity of the galactic plane. We also show results for a 
more conservative case using MWA5000 over three fields 
of view rather than 16. 

For all surveys we assume a redshift range correspond- 
ing to a factor of 3 in frequency and break each observa- 
tion into redshift bins of width Az — 0.1(1 + z). The fac- 
tor of 3 corresponds to the largest frequency bandwidth 
over which a low frequency dipole antenna has suitable 
sensitivity. In our lower redshift case we have assumed 
that dishes will be used. These could potentially have a 
larger bandwidth than the dipole antennae, but we con- 
servatively use the same factor of 3 in frequency. We 
also assume that T sys = T s k y + Tinst where the sky tem- 
perature T sky = 260[(l + z)/9.5] 2 - 55 >2], and that the 
instrumental temperature is 30K in an optimistic case 
and 100K in a pessimistic case, corresponding to what is 
reasonable with current technology and what might be 
expected by the time the hypothetical experiments we 
discuss will be built. 
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B. Cosmological Parameters 

Throughout this work we assume that the true back- 
ground cosmology is that of a flat ACDM Universe 
with density parameters Q\ = 0.7 in dark energy 
Q m h 2 = 0.147 in matter, Q h h 2 = 0.023 in baryons, 
and f2„/i 2 = 0.00054 in neutrinos, and with values of 
A* = 25.0 x 10~ 10 , n s = 0.95, and a = 0.0 for the 
spectral index, amplitude, and running of the primordial 
power spectrum. The chosen value of A s corresponds to 
an rms fluctuation amplitude on a 8/i -1 Mpc scale of 
(Ts = 0.84. We have assumed a neutrino hierarchy with 
one dominant species having a neutrino mass of 0.05eV. 
We assume a fiducial value of r = 0.1 for the optical 
depth to electron scattering during reionization, and a 
helium mass fraction from big bang nucleosynthesis of 
Yuc = 0.24. In a lower redshift example we also leave the 
dark energy equation of state as a free parameter with a 
fiducial value of w = — 1. 



V. RESULTS 

We present the results of our Fisher matrix calcu- 
lations in Tables Hl lVIII We consider several different 
sets of assumptions regarding the effects of the ionizing 
background, which as explained above, may introduce 
a weakly scale dependent bias. In principle this could 
produce degeneracies which would weaken constraints on 
cosmological parameters. First we present the case where 
the scale dependence of the bias is insignificant compared 
with the precision of the power spectrum measurements. 
Then we assume that the effect is large enough to require 
using a model with additional parameters to accurately 
fit the data. We also discuss a very pessimistic scenario 
where the form of the bias is completely unknown. 

We calculate constraints for both optimistic and pes- 
simistic antenna temperature, which are 30K and 100K 
respectively. We also show results for a cosmic variance 
limited version of each survey where the detector noise is 
assumed to make no contribution. This case represents 
the best conditions imaginable for constraining cosmo- 
logical parameters given the finite volume of space that 
is observed. As it turns out, an FFTT is nearly cos- 
mic variance limited at the redshifts of interest here. In 
this case, we only show results for the optimistic antenna 
temperature case. 

We consider a flat ACDM cosmology except for one 
case centered around z — 1.5 where we leave w as a free 
parameter. 

A. Constant Bias 

We begin with the case where the ionizing background 
does not introduce any significant modulation to the 
power spectrum. Thus, the bias, B(k) — b can be taken 
as a constant in each redshift bin. Note that the value 



of the bias is different in each redshift bin, by constant 
we are referring to its dependence on scale. This is the 
most optimistic model one could hope for. Here we as- 
sume that the shot noise, neutral fraction and the bias in 
each redshift bin are unknown. Each is marginalized over 
as a free parameter in the Fisher matrix. We use Eqs. 
([3]) and (0| to determine the fiducial values of the bias. 
The shot noise is determined from the Press-Schechter 
mass function of dark matter halos modified to include 
non-spherical collapse. The fiducial value of the neutral 
fraction is set to Xm — 0.02. The results are presented 
in Tables lIVj 



B. Scale Dependent Bias 

We consider the case where the fluctuations of the 
21cm power spectrum due to the ionizing background 
are large enough to require fitting additional nuisance 
parameters. We do not claim to derive the exact form 
of the modulation, instead we assume an idealized form 
and fit, and then see how cosmological constraints are 
degraded as more flexibility is introduced. As discussed 
above the scale dependent bias, B(k), is the product of 
the average galaxy bias, b and (1 — iC(k)). In [17] , K(k) 
is calculated for the case where all photons travel the 
same distance. Here we consider the more realistic case 
of a constant mean-free-path. In this case, as calculated 
in appendix A, the modulation takes the form -K"(k) = 
ifoa r ctan(/cA m fp)/A m f p fc. We use the following approxi- 
mation (given here in proper distance, but appearing in 
comoving distance in K(k)) for the fiducial values of the 
mean- free-path: A mfp = 85((l + z)/4) 4 Mpc H. As 
before, in our Fisher matrix calculations we marginalize 
over shot noise, neutral fraction and average bias in each 
redshift bin. Additionally, we marginalize over Kq and 
A m fp in each redshift bin. Our results are shown in Tables 

EKE 

To see how constraints are affected as we increase the 
flexibility of the fit, we consider a case with small correc- 
tions to the form used above. Here we assume that the 
modulation is of the form 

K (k) = Kq arctan(fcA m f p )/A m fpfc + c\k + c 2 k 2 + c 3 k 3 ... 

(12) 

where the c„'s are set to fiducial values of zero and 
marginalized over in each redshift bin. Aside from 
marginalizing over these new parameters, the Fisher ma- 
trix calculation is done the same as described above. 

Table [*VTT1 shows the results for FFTT for marginaliz- 
ing over all of the c„'s up to degree n. This illustrates 
how parameter constraints will be degraded as there is 
more flexibility in the fit for K(\t). We find that most 
cosmological parameters are not degenerate with scale 
dependent changes to the bias. However, the spectral in- 
dex of the primordial power spectrum and its running are 
degraded significantly. This is because a change in these 
parameters smoothly modifies the 21cm power spectrum 
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as a function of k, which mimics changes due to scale de- 
pendent bias parameters. A change in other parameters 
alters features in the power spectrum that can be dis- 
tinguished from the scale dependent bias. For example, 
one might expect f2„/i 2 to be strongly degenerate as well, 
since neutrinos change the power spectrum in a smooth 
way. However, with our choice of cosmological parame- 
ters, we hold fl m h 2 constant when taking derivatives of 
f2„/i 2 . This means that changing the energy density of 
neutrinos also changes the energy density of cold dark 
matter. The associated change in the cold dark matter 
component alters features in the power spectrum which 
break the degeneracy with the scale dependent bias. 

Note that while the above fit gives a range of how much 
fitting the scale dependent bias could weaken cosmolog- 
ical constraints, future work will have to determine the 
exact form of the fit which should be used. This will 
require detailed numerical simulations. 

We have also considered the extremely pessimistic case 
where the form of the scale dependent bias is completely 
unknown. It has been found that it may be possible to 
isolate the part of the power spectrum containing only 
cosmological information by exploiting the angular de- 
pendence of the power spectrum. We use the PESS ion- 
ization power spectrum model described in [l(| to calcu- 
late Fisher matrices. This isolates the term in the 21cm 
power spectrum which does not depend on the bias. We 
find that even FFTT cannot constrain cosmological ef- 
fectively and does not improve constraints from Planck 
significantly. This is probably overly pessimistic since 
the changes to the bias from the ionizing background are 
expected to be small, but it shows the importance of un- 
derstanding the form of the bias. 



C. Other Factors 

There are several other factors that affect the con- 
straints attainable with post-reionization 21cm surveys. 
We investigate how these factors change our constraints 
in the case of a constant bias. In Fig. [2J we have plotted 
the dependence of the la uncertainty on the bandwidth 
over which foregrounds are removed. Larger values for 
the bandwidth allow measurement of smaller k values 
which are otherwise not accessible. Overall, the sensitiv- 
ity to the bandwidth variation is relatively modest. 

In Fig. [3l we show how the constraints vary for dif- 
ferent choices of fc max , the maximum wavenumber being 
used. We show how the parameter uncertainties change 
if we impose the cutoff at different levels of discrepancy 
between the non-linear and linear power spectra. As 
&max is lowered errors on parameters are substantially 
increased. For a discrepancy of 0.1 and 0.05, /c max ranges 
from 0.68-0. 16MPC" 1 and 0.28-0.12Mpc _1 respectively, 
see Fig. [TJ For a discrepancy of 0.01 fc max s» 0.06Mpc _1 
for each redshift bin. 

In Fig. 21 we show how the constraints depend on 
M m j n , the minimum halo mass above which significant 



1.5 



1.4 - 




0.9 1 1 1 1 1 

10 15 20 25 30 

Bandwidth[MHz] 

FIG. 2: Sensitivity of the relative la error bar to the band- 
width over which foregrounds can be removed for FFTT. Er- 
rors have been normalized such that they equal unity for the 
assumptions in Table [I] The dashed line gives the error for 
Q.A, the dotted line for f!„/i 2 , the solid line for n B and the 
dot-dash line for a. 




Discrepancy 



FIG. 3: Sensitivity of the relative la error bar to the maximal 
wavenumber fc max for FFTT. Errors have been normalized 
such that they equal one for the assumptions in Table U The 
dashed line gives the error for Q,a, the dotted line for fi„/i 2 , 
the solid line for n s and the dot-dash line for a. 



amounts of neutral hydrogen may assemble. A high value 
of M m i n increases the shot noise and degrades the con- 
straints. We find that the shot noise only affects the 
results significantly above ~ 10 12 M Q , which is unrealis- 
tic since hydro gen is found in lower mass galaxies in the 
local Universe I54II. 
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TABLE I: Constant Bias: la Errors on cosmological parameters with surveys centered at z = 3.5. OPT refers to an antenna 
temperature of 30K, PESS 100K, and CVL to the cosmic variance limited case. OPT and PESS correspond to what may be 
possible by the time these experiments will be built and what is achievable today. 







Q A 


r2 m /i 2 


C2 h h 2 




A 2 x 10 10 




n u h 2 




*Hc 


Fiducial Values 




0.7 


0.147 


0.023 


0.95 


25.0 


0.0 


0.00054 


0.10 


0.24 


MWASk-Hcmisphcro 


PESS 
OPT 
CVL 


0.0041 
0.0024 
0.00059 


0.0098 
0.0041 
0.00025 


0.0024 

0.001 

0.000074 


0.03 

0.014 

0.003 




0.012 

0.0068 

0.0013 


0.0018 

0.00095 

0.0003 






MWA5k-3 Fields 


PESS 
OPT 
CVL 


0.0095 
0.0055 
0.0014 


0.023 

0.0096 

0.00059 


0.0056 
0.0023 
0.00017 


0.07 

0.034 

0.007 




0.029 
0.016 
0.003 


0.0042 
0.0022 
0.00071 






FFTT 




0.00076 


0.00032 


0.000087 


0.0034 




0.0016 


0.00032 






Planck 

+ MWA5k-Hemisphere 
+ MWA5k-3 Fields 
+FFTT 




0.038 
0.0023 
0.0053 
0.00076 


0.0041 
0.00044 
0.00074 
0.00022 


0.00024 
0.00013 
0.00017 
0.000076 


0.0094 
0.0037 
0.0056 
0.0027 


0.25 
0.21 
0.22 
0.2 


0.007 
0.0028 
0.004 
0.0013 


0.0039 
0.00041 
0.00077 
0.00026 


0.0046 
0.004 
0.0041 
0.004 


0.014 
0.0057 
0.0082 
0.0043 


a When not otherwise indicated the "OPT" 
used. 


value of T sys has been 














TABLE II: Constant Bias: la Errors on cosmological parameters with surveys centered at z — 1.5. OPT refers to an antenna 
temperature of 30K, PESS 100K, and CVL to the cosmic variance limited case. OPT and PESS correspond to what may be 
possible by the time these experiments will be built and what is achievable today. 






f2 A 




n h h 2 




A 2 X 10 10 








*He 


Fiducial Values 




0.7 


0.147 


0.023 


0.95 


25.0 


0.0 


0.00054 


0.10 


0.24 


MWA5k-Hcmisphcrc 


PESS 
OPT 
CVL 


0.0033 
0.0023 
0.0018 


0.016 
0.011 
0.0085 


0.004 

0.0027 

0.0021 


0.051 
0.036 
0.029 




0.019 
0.014 
0.011 


0.0027 

0.002 

0.0016 






MWA5k-3 Fields 


PESS 
OPT 
CVL 


0.0077 
0.0054 
0.0042 


0.038 
0.026 
0.02 


0.0093 
0.0063 
0.0048 


0.12 

0.084 

0.068 




0.044 
0.032 
0.026 


0.0063 
0.0046 
0.0038 






Planck 

+ MWA5k-Hemisphere 
+ MWA5k-3 Fields 




0.038 

0.0021 

0.0047 


0.0041 

0.00042 

0.00068 


0.00024 
0.00014 
0.00018 


0.0094 
0.0042 
0.0064 


0.25 
0.22 
0.22 


0.007 

0.0033 

0.0048 


0.0039 

0.00041 

0.00075 


0.0046 
0.0041 
0.0041 


0.014 

0.0065 
0.0097 



isolate the term which only depends on cosmology can- 
not effectively constrain cosmological parameters. For 
this reason, it is crucial that the modulation of the 21cm 
power spectrum from the ionizing background is modeled 
reliably. 

Our results show approximately how constraints would 
be degraded if a fit were to be necessary. Table I VIII il- 
lustrates how more freedom in the fit affects cosmologi- 
cal constraints. We find that the only parameters which 
are strongly degenerate with the new fitting parameters 
are n s and a, the spectral index and the running of the 
primordial power spectrum. The la errors on these pa- 
rameters could increase by an order of magnitude. The 
other la cosmological parameter errors are increased by 
roughly 20% or less. 

For comparison, the constraints from 21cm surveys 
that focus on the EoR vary considerably under differ- 
ent ionization power spectrum model assumptions. The 
robustness of our post-reionization constraints is mainly 
due to the fact that we do not have to consider the ioniza- 
tion power spectrum, only small scale dependent changes 
to the bias. 

Our results also depend strongly on the maximum 
wavenumber of the power spectrum which can be used 
to reliably measure cosmological parameters, k max . It 
will be desirable for theoretical predictions of the matter 
power spectrum, presumably from N-body simulations, 
to be computed to a fraction of a percent for k max val- 
ues corresponding to a 5% descrepancy in the linear and 
non-linear power spectra shown in Fig. [1] If the matter 



used. 




M . 

mm 



FIG. 4: Dependence of constraints on M m in (in units of solar 
mass) for FFTT. Errors have been normalized such that they 
equal one for the assumptions in table [I] The dashed line 
gives the error for Qa, the dotted line for fi„/i 2 , the solid line 
for n s and the dot-dash line for a. 



VI. DISCUSSION 



The above results are very promising. We find that 
using the angular dependence of the power spectrum to 
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TABLE III: Ratio of la Errors for surveys combined with Planck to Errors from Planck alone. 



f> A S"2 m h 2 Q b h 2 Tig A 2 a £l v h? t V Hc 



MWASk 


Hei 


liaphcrc 




= 3 




0.061 





11 


0.54 


0.39 


0.84 


0.4 


0.11 


0.87 


0.41 


MWA5k 


3 F 


elds 




= 3 




0.14 





18 


0.71 


0.6 


0.88 


0.57 


0.2 


0.89 


0.59 


FFTT 








= 3 


5 


0.02 





054 


0.32 


0.29 


0.8 


0.19 


0.067 


0.87 


0.31 


MWASk 


Hei 


liaphcrc 




= 1 




0.055 





1 


0.58 


0.45 


0.88 


0.47 


0.11 


0.89 


0.46 


MWA5k 


3 F 


elds 




= 1 




0.12 





17 


0.75 


0.68 


0.88 


0.69 


0.19 


0.89 


0.69 



"The "OPT" value of T sya (30K) and the constant bias model 
have been used. 



TABLE IV: Constant Bias: la Errors on cosmological parameters with surveys centered at z — 1.5 with w as a free parameter. 



OPT refers to an antenna temperature of 30K, PESS 100K, and CVL to the cosmic variance limited case. OPT and PESS 
correspond to what may be possible by the time these experiments will be built and what is achievable today. 






"A 


f! m fi 2 


S2 b h 2 




A 2 X 10 10 




n„h 2 








Fiducial Values 




0.7 


0.147 


0.023 


0.95 


25.0 


0.0 


0.00054 


-1.0 


0.10 


0.24 


MWA5k- Hemisphere 


PESS 
OPT 
CVL 


0.0052 
0.0047 
0.0046 


0.017 
0.012 
0.0088 


0.004 

0.0027 

0.0021 


0.051 
0.036 
0.029 




0.019 
0.014 
0.012 


0.003 

0.0022 

0.0018 


0.047 
0.039 
0.036 






MWA5k-3 Ficlda 


PESS 
OPT 
CVL 


0.012 
0.011 
0.011 


0.039 
0.027 
0.02 


0.0094 
0.0064 
0.0048 


0.12 

0.084 

0.068 




0.045 
0.033 
0.027 


0.0069 
0.0051 
0.0043 


0.11 

0.092 

0.083 






Planck 

4- M WA5k-Hemisphere 
4-MWA5k-3 Fields 




0.096 

0.0038 

0.0082 


0.0061 

0.00056 

0.00097 


0.00024 
0.00014 
0.00018 


0.0094 
0.0042 
0.0064 


0.27 
0.22 
0.22 


0.0071 
0.0033 
0.0048 


0.0059 

0.00049 

0.0009 


0.16 

0.023 
0.045 


0.0051 
0.0041 
0.0041 


0.015 

0.0065 

0.0097 



"When not otherwise indicated the "OPT" value of T sys has been 



used. 



TABLE V: Fit to Scale Dependent Bias: la Errors on cosmological parameters with surveys centered at z — 3.5. OPT refers 
to an antenna temperature of 30K, PESS 100K, and CVL to the cosmic variance limited case. OPT and PESS correspond to 
what may be possible by the time these experiments will be built and what is achievable today. 







n A 


n m h 2 


n h h 2 




A 2 x 10 10 




£l u h 2 




*Hc 


Fiducial Values 




0.7 


0.147 


0.023 


0.95 


25.0 


0.0 


0.00054 


0.10 


0.24 


MWASk-Hcmisphcrc 


PESS 


0.0043 


0.011 


0.0027 


0.047 




0.028 


0.002 








OPT 


0.0025 


0.0046 


0.0011 


0.028 




0.017 


0.001 








CVL 


0.00061 


0.00026 


0.000077 


0.007 




0.0031 


0.00033 






MWA5k-3 Fields 


PESS 


0.01 


0.026 


0.0062 


0.11 




0.066 


0.0047 








OPT 


0.0058 


0.011 


0.0026 


0.065 




0.04 


0.0024 








CVL 


0.0014 


0.0006 


0.00018 


0.016 




0.0072 


0.00078 






FFTT 




0.00079 


0.00033 


0.000089 


0.0085 




0.0041 


0.00035 






Planck 




0.038 


0.0041 


0.00024 


0.0094 


0.25 


0.007 


0.0039 


0.0046 


0.014 


+ MWA5k-Hemisphere 




0.0024 


0.00052 


0.00013 


0.0042 


0.22 


0.0032 


0.00051 


0.0041 


0.007 


+ MWA5k-3 Fields 




0.0055 


0.00078 


0.00017 


0.0063 


0.22 


0.0048 


0.00094 


0.0041 


0.01 


+ FFTT 




0.00079 


0.00025 


0.000078 


0.0034 


0.21 


0.0017 


0.00029 


0.004 


0.0052 




a When not otherwise indicated the "OPT" 


value of T s 


ys has been 















used. 



TABLE VI: Fit to Scale Dependent Bias: la Errors on cosmological parameters with surveys centered at z — 1.5. OPT refers 
to an antenna temperature of 30K, PESS 100K, and CVL to the cosmic variance limited case. OPT and PESS correspond to 
what may be possible by the time these experiments will be built and what is achievable today. 







f2 A 


fl m h 2 


n h n 2 




A 2 x 10 10 




£l u h 2 




*He 


Fiducial Values 




0.7 


0.147 


0.023 


0.95 


25.0 


0.0 


0.00054 


0.10 


0.24 


MWASk-Hcmisphcrc 


PESS 
OPT 
CVL 


0.0034 
0.0024 
0.0019 


0.019 
0.013 
0.0095 


0.0047 
0.0031 
0.0023 


0.067 
0.048 
0.039 




0.034 
0.025 
0.02 


0.0032 
0.0023 
0.0019 






MWA5k-3 Fields 


PESS 
OPT 
CVL 


0.008 

0.0056 

0.0043 


0.045 

0.03 

0.022 


0.011 

0.0072 

0.0053 


0.16 
0.11 
0.09 




0.08 

0.058 

0.046 


0.0076 
0.0054 
0.0044 






Planck 

+ MWA5k-Hemisphere 
+ MWA5k-3 Fields 




0.038 

0.0021 

0.0048 


0.0041 

0.00049 

0.00071 


0.00024 
0.00014 
0.00018 


0.0094 
0.0045 
0.0069 


0.25 
0.22 
0.22 


0.007 

0.0037 

0.0052 


0.0039 

0.00053 

0.00092 


0.0046 
0.0041 
0.0042 


0.014 

0.0073 

0.011 


a When not otherwise indicated the "OPT" 


value of T sy 


3 has been 















used. 



TABLE VII: Polynomial Corrections to Bias: la Errors on cosmological parameters for FFTT centered around z = 3.5. Effects 
of fitting small polynomial corrections to the scale dependent bias. Here n is the degree of the polynomial marginalized over. 
The case with n=0 corresponds to the value from table IVl 



constant 

n=0 

n=l 

n=2 

n=3 

n=4 



0.7 

0.00076 
0.00079 
0.00084 
0.00092 
0.00093 
0.00093 



0.147 

0.00032 

0.00033 

0.00033 

0.00035 

0.00037 

0.00038 



0.023 

0.000087 

0.000089 

0.00009 

0.000093 

0.000097 

0.0001 



0.95 

0.0034 

0.0085 

0.011 

0.016 

0.017 

0.017 



0.0 

0.0016 

0.0041 

0.0071 

0.014 

0.015 

0.015 



0.00054 
0.00032 
0.00035 
0.00036 
0.00036 
0.00037 
0.00037 
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power spectrum can not be computed to this accuracy, 
the value of fc max may be the limiting factor for measur- 
ing cosmological parameters with experiments similar to 
those described in this paper. 

The noise in the interferometers is lower at lower red- 
shift due to the redshift scaling of the galactic syn- 
chrotron emission. This advantage is offset by the fact 
that at higher redshifts the neutral fraction is higher and 
the observed power spectrum stays linear up to a higher 
fc max value. The FFTT and the MWA5000 at z = 1.5 
are nearly cosmic variance limited. These experiments 
would be close to capturing all the information available 
from the power spectrum available to us on linear scales 
at the relevant redshifts. The MWA5000 with optimistic 
antenna temperature centered at z = 3.5 gives errors 
roughly five times larger than the cosmic variance lim- 
ited case. 

Current cosmological information sets the upper limit 
on neutrino mass to ~ 0.62 eV [H, H^. For FFTT 
with a constant bias, we find a la error corresponding 
to 0.030eV. As expected [l5|, this is clearly much bet- 
ter than existing constraints and is sufficiently small to 
probe the neutrino mass differences implied by neutrino 
oscillations. 

It is also interesting to compare how well these surveys 
constrain cosmology with respect to other future probes. 
The FFTT survey will perform much better than the 
Planck satellite for almost all of the parameters under 
consideration. When combined with Planck, essentially 
all of the constraints are improved. In Table IIIII we ex- 
plicitly show how the constraint on each parameter from 
Planck would improve if combined with post-reionization 
21cm surveys. The most improved parameters are f?A, 
fl m h 2 , and Q„/i 2 , which for FFTT are 2.0%, 5.4%, and 
6.7% of the errors from Planck alone. 

Since most of the remaining neutral hydrogen after the 
EoR is found in DLAs hosted by galaxies, we expect our 
results to be similar to galaxy redshift surveys over sim- 
ilar redshift ranges. In [1 11 ] an all sky galaxy survey cen- 
tered around z = 1, termed G3, is considered. We find 
the cosmological constraints they calculate to be very 
similar to those in our MWA5000 case with a constant 
bias centered around z = 1.5. They also consider a Ly- 
man break survey centered around z = 3, which they 
term G2. This survey contains roughly 1% of the vol- 
ume included in our MWA5000-hemisphere and FFTT 
cases. Since the error on the power spectrum scales as 
the square root of the volume of a survey, this should 
cause parameter errors in G2 to be roughly an order of 
magnitude larger. We find that this falls roughly into the 
range we calculate above for the constant bias MWA5000 
case centered around z = 3.5. Note that if the degen- 
eracies due to the UV background discussed above are 
present, galaxy surveys may do significantly better than 
21cm at measuring the n s and a, the spectral index of 
the primordial power spectrum and its running. 

In the case where we allow w to vary as a free pa- 
rameter we find the constraints on most parameters es- 



sentially unchanged. The main exception is Q\, where 
we find that our constraints are increased by roughly a 
factor of two. The constraints achievable for w are inter- 
esting, about 3.9% if the whole hemisphere is observed, 
and 9.2% with only three fields of view. 



VII. CONCLUSIONS 

A 21cm survey after reionization provides a promising 
probe of cosmology and fundamental physics [l5j . We 
find the simple form of the 21cm power spectrum after 
reionization provides less uncertainty in forecasting con- 
straints than 21cm during the EoR. However, we find that 
errors on estimates of the spectral index and the running 
of the primordial power spectrum could be increased by 
an order of magnitude if the ionizing background intro- 
duces scale dependent effects into the post reionization 
21cm power spectrum. 

Aside from its greater simplicity, the measurement sen- 
sitivity after reionization is improved over EoR surveys 
due to lower foreground brightness temperature, a bias 
that is greater than unity, and a larger growth factor of 
density perturbations. These advantages are somewhat 
offset by the fact that there is less neutral hydrogen after 
reionization. 

The constraints we derive are competitive with those 
from the EoR and for the FFTT are significantly better 
than next generation CMB experiments on their own. 
We find errors on SIa, ^m^ 2 , and Vl u h 2 to be improved 
the most, each lowering the errors from Planck alone by 
at least an order of magnitude for both FFTT and the 
survey similar in design to MWA5000. 



APPENDIX: MODULATION OF 21CM PS FROM 
IONIZING BACKGROUND 

We derive the form of the modulation of the 21cm 
power spectrum due to the ionizing background in the 
case of a constant mean- free- path. The only difference in 
the derivation presented in [17| is the form of the fluctu- 
ations of the ionizing flux <5j(x). We convolve the real- 
space density field with a filter function to account for the 
effects of a finite mean-free-path and the inverse squared 
dependence of ionizing flux. The fluctuation in flux is 
given by 



<5j(x) oc / dx'G(x,x')<5(x'). (A.l) 

We assume a mean-free-path of ionizing photons which 
is independent of position such that, 

G(X)XOoc !fE(z^i|W ) {A2) 
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where A m f p is the mean free path of ionizing photons. 
The convolution theorem gives 

6j(k) oc g(k)S k , (A.3) 

where g(k) and Sk are the Fourier transforms of the filter 
function and the fluctuation in ionizing flux. We find 
that 

g(k) oc arctan(fcA m f p )/fc. (A. 4) 

From the derivation presented in [TtJ , this implies that 
the power spectrum is given by Eqn. [2] with 



K(k) = A'oarctan(fcA m f p )/A m fpfc, (A. 5) 

where Kq is a constant of the order of 0.01 or smaller 
0- 
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